Two-electron bound states in continuum in quantum dots 
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Bound state in continuum (BIG) might appear in open quantum dots for variation of the dot's 
shape. By means of the equations of motion of Green functions we investigate effect of strong 
intradot Goulomb interactions on that phenomenon in the framework of impurity Anderson modeL 
Equation that imaginary part of poles of the Green function equals to zero gives condition for BIGs. 
As a result we show that Goulomb interactions replicate the single-electron BIGs into two-electron 
ones. 

PACS numbers: 73.23.Hk, 72.10.-d, 71.27.-l-a, 05.60.Gg 

In 1929, von Neumann and Wigner [l| predicted the existence of discrete solutions of the single-particle Schrodinger 
equation embedded in the continuum of positive energy states. Their analysis examined by Stillinger and Herrick [3| 
long time was regarded as mathematical curiosity because of certain spatially oscillating central symmetric potentials. 
More recently in 1973 Herrik ^] and Stillinger [1] predicted BICs in semiconductor heterostructure superlattices which 
were observed by Capasso et al as the very narrow absorption peak Q . 

In the framework of the Feshbach's theory of resonances Friedrich and Wintgen Q have shown that BIC occurs 
due to the interference of resonances. If two resonances pass each other as a function of a continuous parameter, then 
for a given value of the parameter one resonance has exactly vanishing width. Later this result was reproduced in 
application to different physical systems in the two-level approximation 0, H, H, [lOj [iHj [Hi • Straight waveguide with 
an attractive, finite size impurity presents an example of realistic structure in which Kim et al presented the numerical 
evidence for the BIC for the variation of the impurity size [T3|. Further, calculations in microwave and semiconductor 
open structures showed that the resonance width also can turn to zero for variation of angle of bent waveguide , 
shape of quantum dot (or resonator) [l^i [iBl , or magnetic field f]3. Recently it was rigorously shown that the zero 
resonance width is the necessary and sufficient condition for BIC (l6l [TtI . This condition means that a coupling of 
the resonance state with continuum equals zero to convert the state into BIC [l^ [l^ . 

That very restricted list of references shows that BICs might occur in different open quantum systems including, for 
example, laser induced continuum structures in atoms p^ . in the molecular system j20| . However in application to 
open quantum dots (QD) BICs were studied in the single electron approximation whereas the Coulomb interactions 
between electrons might be very important for robustness of the BIC. In the present work we consider effect of local 
intradot Coulomb interactions in QD onto BICs in the framework of the two-level impurity Anderson model plj that 
is one of the most important theoretical models for a study of strong correlations in condensed matter physics. 

We consider quantum dot (QD) coupled to leads (left and right) which support one propagating mode (the case of 
two continuums) with the following total Hamiltonian 

^ Hc + Hd + V. (1) 

C=L,R. 

The leads, left (L) and right (R) in ([T]) are presented as the non interacting electron gas 

Hc = Yl e(k)c+^cckac. C^L,R. (2) 

A continual spectrum e(fc) defines the propagating band of leads. The Hamiltonian of many level QD is that of the 
impurity Anderson model [2l| . 
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Here a+g. is the creation operator of an electron on the m-th level of QD, Um takes into account the Hubbard repulsion 
at the level to, and rima = o-ma^rna- The interaction 

Vm{k){clca,na + h.c) (4) 

kamC 

describes couplings between the leads and QD where c^^^ is the creation operator of an electron in the leads C. 
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In order to calculate transport properties of the QD we use a technique of the equations of motion for retarded and 
advanced Green functions which successfully used to consider the Fano and Kondo resonances in the Anderson model 
[H, m, [IE, HE US] ■ Following Laxroix [28| we use a Hartree-Fock approximation in the wires {{ckaC o-naC'na \ a'^a' ) ) ~ 
(i^na) {{ckcrcWma')) ■ approximation is justified for weak couplings compared to the Coulomb interactions: V,n ^ 
Um- As a result we obtain the following equation 

G-\E) = GQ\,iE) + tr (5) 

for the Green functions Gma,na'{E) = {{0'ma\ci^cr'))^^ ^^^'^ form of the Dyson equation '24*1 . Here Gqd{E) is the 
Green function of the isolated QD 

GQD,mm' ,a,cr'{E) — GQD,mcr{E)S 

mm' ^a.<7' 

GQD,raAE) = + E^^^- (6) 

which are exact for the isolated QD. For the simplicity we take wide band wires and approximate the self-energy as 

MM 



E 



k 



E-e{k)^ + iO 



= -inV„iVnPc{E) = -ivT^nT^ (7) 



where pc{E) is the density of states of the left and right wires. The average values of the occupation numbers 
fa^o-} = {O'maO'ina) which enter the expressions for the Green functions are calculated self-consistently via the formulas 

{nma) J d-EImGrna.rnaiE). (8) 

The form of the self-energy ^ and the QD Green function ([6]) allows to proceed to the case of free electrons with 
Um — 0. In this case BIG appears if QD acquires accidental degeneracy ei = £2 Therefore in the vicinity of 

the degeneracy point e = £2 — ei = we can restrict ourselves to the two-level approximation for QD 9]. Then the 
occupation numbers ((S]) are given by four poles of the Green fimction At zero temperature, the transmission 
amplitude can be expressed in term of the Green function 

r-rG(s)r+, r = (ri,r2). (9) 

The results of numerical self-consistent calculation of the transmission ^ are presented in Fig. [T] Fig. [T] (a) shows 
the case of zero Hubbard repulsion Umn = (no electron correlations) in which QD is given only by the single-electron 
energy levels. As shown in BIG occurs here at the point of degeneracy of electron states in QD for e = 0. At this 

point the S-matrix becomes singular because the transmission zero crosses the unit transmission [12] as shown in Fig. 
[T]where the unit transmission follows the energy levels. As the Hubbard repulsion is included, QD is given not only by 
single-electron states but also by two-electron states as shown in Fig. [l](b) by solid lines. As a result we obtain that 
the number of degenerated points becomes four as seen from Fig. [1] (b). One can see that lines of zero transmission 
cross the lines of maximal unit transmission at these points. Therefore, one can expect the BICs at four points of 
degeneracy of the QD. In order to show that we present in Fig. [2] the resonance widths of the energy levels defined as 
Fa = —2Im{z\), A = 1, 2, 3, 4, where z\{E, e) are the poles of the Green function or zeros of the right hand expression 
in the Dyson equation ([S]). The points at which F^ = define BIGs [H, [l3- One can see that these points coincide 
with the points of degeneracy of the QD given by equations ed = 0, — Ui/2, = —U2/2, ec4 = {U2 — Ui)/2. 
The corresponding energies of BIGs are 0, 0.1, — 0.15, 0.05. The first BIG is pure single electron localized state 
superposed of two single electron states with m = 1,2 [l2|. However the next two BIGs are supespositions of the 
single electron states and two-electron ones. At last, for the last case (ec4 = {U2 — Ui)/20 the BIG is superposed of 
the two-electron states in QD. Although specific values of F™ has no importance for BIG's points defined by crossings 
of the energ;^levels of QD, they are important for appearance of the Dicke superradiant state which accumulates the 
total width [9] as seen from Fig. [51 

Since the resonance width turns to zero with approaching the BIG point, we expect singular behavior of occupation 
numbers ^ at the energy of BIG. In fact. Fig. [5] (a, b, c) demonstrate this effect. Let us consider the first BIG at 
e = with discrete energy E = Q dX which the single-electron energies in QD are crossing as shown in Fig. [T] (b). 
One can see from Fig. [3] (a) that at the energy E=0 both energy levels are sharply and simultaneously populated till 
one half. The next resonances with finite widths correspond to the two-electron energies of QD that are populated 
smoothly at the Hubbard repulsive energies Ui = 0.2 and U2 = 0.3 by usual scenario as seen from Fig. [3] (a). 
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FIG. 1: (a) The transmission In |r| of QD versus energy of incident electron and energy splitting e for the case of zero Ifubbard 
repulsion U\ = U2 = 0. (b) The case of strong Hubbard repulsion Ui — 0.2, U2 ~ 0.3, Pi = r2 = 0.05. The single-electron 
and two-electron energy levels in closed QD are shown by thin lines. Black regions correspond to those where the transmission 
close to zero while the white regions do to the maximal transmission. 
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FIG. 2: The resonance widths defined as —2Im[z\{E, e)], A = 1, 2, 3, 4 versus e for E = where z\ are the poles of the Green 
function ((5]). 

The next BICs happen for the single-electron state crosses the two-electron state at points e = —0.15 and e = 0.1. 
Because of similarity of these BIC points we have presented here only the first case as shown in Fig. [3] (b). The 
BIC's discrete energy for that case equals to E — 0.15 (Fig. [1]). Again we see that for approaching to this energy the 
BIC populates sharply. However the populations of the single-electron level and two-electron one are well separated 
because of the Hubbard repulsion of the two-electron state. The last figure Fig. [3] (c) refers to the crossing of two- 
electron energies at the point e — —0.05. As seen from Fig. [T] (a) the two-electron BIC has energy E — 0.25. As a 
result for approaching this energy we observe sharp population of this state similar to the case (a). 

Are BICs critical to energy level crossing? Similar to [1, [l^] we lift the degeneracy in the QD by transitions 
between levels, adding a hopping term between the single-electron states into the Hamiltonian of the two-level QD, 
Hd — » Hd — va^^a^ij — va^^aia which evolves the picture of energy crossing into a picture with an avoided crossing. 
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FIG. 3: The electron populations as dependent on the energy of incident electron defined by ((Sjl for the parameters of the 
system given in Fig. [T] (a) e = 0, (b) e = —0.15, and (c) e = —0.05. The transmission \T\ is shown by thin green line. 
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FIG. 4: (a) The transmission ln(— ln(l — \T\)) of QD versus energy of incident electron and energy splitting e in the avoiding 
crossing scenario v — 0.05, Ui — 0.2, U2 ~ 0.3, Fi = r2 = 0.05. Black regions correspond to those where the transmission close 
to zero while the white ones do to those where the transmission is near unit. Thereby the white regions follow the single-electron 
and two-electron energy levels in closed QD. Positions of BICs are shown by open circles, (b) The resonance widths dependent 
on e. 

Fig. m (a) shows the transmission of QD in which the energy levels repel each other because of the hopping between 
QD levels. In order to show clearly the zero and unit transmission we present the transmission in the double log scale 
ln(— ln(l — |r|)). One can see the avoided level crossings shown by white lines with T = 1. BICs shown by open 
circles are located at those points where the unit transmission T = 1 (white lines) crosses the zero one G = (black 
lines) similar to the case of non interacting electrons Fig. [l](b) shows that the resonance widths turn to zero at 
four critical values of e. 

The Hubbard repulsion presented in the Anderson impurity model ([3]) is not the only way to account the Coulomb 
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interactions. The last also induce the inter-level couplings in the form '^mn'^a-cr' Umnnmanna' ■ Therefore, in the 
two- level approach a new Coulomb constant U12 appears. The equations of motion for the Green functions in the QD 
become tedious but still complete to give the following Green function 

, . ^ {E + Em){E + Em - Urn - Ul2{nma)){E + Em - Ul2 - U,n{nrn)) 

^""'""^ ' {E + Em - Um{l - {Uma)) - Ui2{nma)){E + Em - t/l2(l - {Uma)) - Um{nma)){E + Em - Um - Ui^iflma))) ' 

(10) 

where m — l,2,m — 2,1, Em — Te- A substitution of (fTU)) into Eqs. (O, ([H]) and ^ allows to calculate the 
transmission of the QD presented in Fig. [5] (a) . Each crossing of energy levels shown by solid lines gives rise to BICs 




FIG. 5: The same as in Fig. |4]but with inclusion of the inter-level Coulomb interactions U\2 = 0.1. (a) The transmission In |T 
of QD versus energy of incident electron and energy splitting e in the avoiding crossing scenario v — 0.05, Ui — U2 = 0.2, Fi — 
F2 — 0.1. (b) The resonance widths dependent on e. 



as shown in Fig. [5] (b). One can see that for e = there are simultaneously four crossings. As a result at this points 
four resonance width turn to zero as shown in Fig. [5](b). Corresponding at the points e = ±0.05 we obtain three BICs 
and so on. If to compare all Figures with transmission probability through the QD one can see that the Coulomb 
interactions in QD replicate the transmission zeros which are between neighboring resonances. If the resonances are 
crossing by an effect of gate voltage we observe BICs at each crossing points Q as Figs [21 [H (b) and [5] (b) show. 
Acknowledgements AF thanks Igor Sandalov for helpful discussions. 
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